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ABSTRACT 

Magnetic-field generation by a relativistic ion beam propagating through an electron-ion plasma 
along a homogeneous magnetic field is investigated with 2.5D high-resolution particle-in-cell (PIC) 
simulations. The studies test predictions of a strong amplification of short-wavelength modes of mag- 
netic turbulence upstream of nonrelativistic and relativistic parallel shocks associated with supernova 
remnants, jets of active galactic nuclei, and gamma-ray bursts. We find good agreement in the prop- 
erties of the turbulence observed in our simulations compared with the dispersion relation calculated 
for linear waves with arbitrary orientation of k. Depending on the parameters, the backreaction on 
the ion beam leads to filamentation of the ambient plasma and the beam, which in turn influences 
the properties of the magnetic turbulence. For mildly- and ultra-rclativistic beams, the instability 
saturates at field amplitudes a few times larger than the homogeneous magnetic field strength. This 
result matches our recent studies of nonrelativistically drifting, hot cosmic-ray particles upstream of 
supernova-remnant shocks which indicated only a moderate magnetic-field amplification by nonrcs- 
onant instabilities. We also demonstrate that the aperiodic turbulence generated by the beam can 
provide efficient particle scattering with a rate compatible with Bohm diffusion. Representing the 
ion beam as a constant external current, i.e. excluding a backreaction of the magnetic turbulence 
on the beam, we observe non-resonant parallel modes with wavelength and growth rate as predicted 
by analytic calculations. In this unrealistic setup the magnetic field is amplified to amplitudes far 
exceeding the homogeneous field, as observed in recent MHD and PIC simulations. 
Subject headings: acceleration of particles, cosmic rays, gamma rays: bursts, methods: numerical, 
shock waves, turbulence 



1. INTRODUCTION 

Collisionless shocks are the acceleration sites of ener- 
getic particles responsible for high-energy emission of as- 
trophysical objects and contributing to the flux of cosmic 
rays detected at Earth. Nonthermal particle populations 
at nonrelativistic shocks of supernova remnants (SNRs) 
are believed to be generated by diffusive shock accel- 
eration (DSA). The particle spectra thus produced in- 
deed agree well with those deduced fro m radio-to-X-ra y 
electron synchrotron emission of SNRs (|Revnolddl20"0l . 
The high efficiency of the DSA mechanism together with 
considerations of the global energetics of a supernova 
explosion make the forward shocks of shell-type SNRs 
the prime candidates for the sources of Galactic cos- 
mic rays (CRs). Nonthermal power-law particle spec- 
tra attributed to relativistic electrons are also inferred 
from modeling the electromagnetic emission from astro- 
physical sources harboring relativistic shock waves, such 
as jets of active galactic nuclei (AGN) and gamma-ray 
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bursts (GRBs) (|Meszarosll20"ol . 

Efficient particle acceleration at shocks invariably re- 
quires a continuous excitation of magnetic turbulence in 
the upstream region, which serves as a scattering medium 
to confine the energe tic particles to the shock r egion for 
further acceleration (jMalkov fc Diamond! 1200 II ). Turbu- 
lent magnetic fields of amplitude much larger than the 
homogeneous interstellar field are needed upstream of 
SNR shocks to account for protons of energies up to 
and beyond the "knee" at 10 15 eV in the cosmic-ray 
spectrum. Recent X-ray observations of several young 
SNRs give evidence that indeed highly ampli fied fields 
exist d ownstream of SNR forward shocks (see IRevnoldi 
(2008), and references therein). Downstream of ultrarel- 
ativistic GRB external shocks the magnetic fields must 
also be amplified orders of magnitude beyond shock- 
compression levels to explain GRB afterglow spectra and 
light curves. Even in the preshock medium magnetic 
fields of milligauss strengths a re required to accou nt for 
the observed X-ray afterglows (|Li fc W axman l2006f) . 

A plausible scenario for magnetic-field generation as- 
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sumes that cosmic ray particles accelerated at the shock 
drift as an ensemble relative to the upstream plasma 
and trigger a variety of instabilities that may lead to 
the g r owth of a t urbulent field component (jBell fe Lucekl 
120011 iBel 120041 ). The distribution function of the cos- 
mic rays is shaped by the scattering rate in the self- 
excited field upstream, thus forcing a nonlinear relation- 
ship between the upstream plasma, the energetic par- 
ticles, and small-scale electromagnetic fields. The up- 
stream field would subsequently be advected and com- 
pressed downstream of the shock and possibly further 
amplified by fluid instabilities in the downstream plasma 
(iGiacalone fc Jokipii|[2007HZirakashvili fc Ptuskinl 120081 : 
iCouch et al.l I2008D . While a full modeling of the up- 
stream region is elusive to date, simulations of turbulence 
build-up using prescribed distribution functions for the 
upstream plasma and the cosmic rays can be invaluable 
tools for the study of the saturation processes and lev- 
els, as well as the backreaction of the evolved turbulence 
on the particles. The system in which a population of 
cosmic rays slowly drifts relative to the upstream plasma 
has been studied with MHD simulations, which rep- 
resented cosmic r ays with a constant ex t ernal current 
(lBell [200l. [20051 : iZirakashvili eTaTI [200l IReville et all 
2008) and with first-principles particle-in-cell (PIC) 
simulations assum ing a constant cosmic-ray current 
(Oh ira et a l. 2009) or including the full d ynamics of the 
energetic particles (iNiemiec et afl 120081 : Istroman et al.1 
120091: iRiquelme fc Spitkovskvll2009h. T hev confirmed the 
quasi- linear predictions by [Bell (2004), who showed that 
for the parameters of young SNRs, magnetic turbulence 
would appear in a form of nonresonant, circularly polar- 
ized, and aperiodic transverse waves. Numerical simu- 
lations analyzing the nonlinear evolution of the system 
found that the turbulence growth eventually saturates, 
but the exact saturation levels differ between the ap- 
proaches, the full PIC simulations typically yielding con- 
siderably lower field amplitudes than MHD studies with 
constant CR currents. 

Here we report kinetic (PIC) simulations of the inter- 
action between the far-upstream plasma and a cold di- 
lute relativistic beam of particles streaming along a ho- 
mogeneous background magnetic field. We assume that 
the beam is composed of cosmic-ray ions, and the cur- 
rent and charge carried by the beam is balanced by 
electrons of the background medium. The situation is 
relevant to the upstream region of both relativistic and 
nonrelativistic shocks undergoing efficient particle ac- 
celeration. In a SNR shock environment, it applies to 
the most energetic cosmic rays accelerated at the shock 
which stream far upstream of the free-escape boundary. 
In this case, a predominantly ionic cosmic-ray compo- 
nent results from the character of injection processes at 
nonrelativistic shocks. Upstream of relativistic shocks of 
GRBs and AGN, distributions of particles accelerated in 
a wide energy range are highly anisotropic. This is be- 
cause particles and the shock move close to the speed 
of light, and a deflection of a particle trajectory by an 
angle greater than TJ^ allows the shock to overtake the 
particle. Therefore, in the upstream rest frame the non- 
thermal particles are highly beamed and their transverse 
momenta are a factor of r"^ 1 smaller than the momenta 
along the shock direction. We approximate this sit- 



uation by assuming that the CR beam is cold. This 
assumption better holds for the freshly accelerated par- 
ticles in the far-upstream region, whose transverse mo- 
menta are much smaller than TJ^ 1 times the parallel mo- 
mentum; the highest-energy particles escape from the 
precursor and therefore have an anisotropic distribution 
in the shock rest frame, which is further enhanced by 
the shock curvature. Some instabilities, e.g. filamenta- 
tion, depend sensitively on the transverse temperature 
of the CR beam, and care must be exercised in extrap- 
olating our simulation results to situations in which the 
CR-beam properties are somewhere between those of the 
cold beam studied h ere and the very hot , but slow beam 
inves t igated earlier dNiemiec et al.l l2008t iStroman et al.1 
2009; Riquel me fc Sr)itkovskvll2009n . However, PIC 
simulations of relativistic shocks in electron-ion plasmas 
suggest that filamentation indeed occurs only far up- 
stream of the shock (Spitkovsky 2008a; 2008b; see also 
Medvedev & Zakutnyaya 2009), and it is generated by a 
warm (p± < p||/r s zj) ion beam. We may therefore ex- 
pect that our assumption of a cold CR beam remains a 
valid approximation for systems with warm CR beams. 
Furthermore, our setup applies to the cosmic-ray ions 
whose energies are larger than the upper limit on the 
energy of electrons accelerated by the s hock, which is 
impo sed by radiative energy losses (e.g. iLi fc Waxmanl 
2006). The highly energetic ions will thus reach farther 
upstream than the CR electrons and the return current 
will be provided by the ambient electrons. 

Note that the applicability of the system under study 
to relativistic astrophysical sources relies on the abil- 
ity of relativistic shocks to accelerate particles to very 
high energies. Although the first-order Fermi process 
at such shocks is widely considered to be the source 
of cosmic ra ys, recent studies in the test particle ap- 
proximation (|Niemiec fc Ostrows ki 2006; Nie miec et all 
Ml: IE emoine et a l. 2006) and using PIC simulations 
( Sironi fc Spitkovskvl l2009f ) show that this mechanism 
can operate only in quasi-parallel or weakly magnetized 
shocks. If the GRB or AGN outflows are strongly mag- 
netized /quasi-perpendicular, some other processes must 
be responsible for particle acceleration (e.g., magnetic 
reconnection), and our results do not apply. 

It is known from studying non-relativistic beams in 
interplanetary space that a competition arises between 
resonant and nonresonant mo des, which exert differ - 
ent backreactions on the beam (Winskc fc Lerovl [l984). 
For the case of a monoenergetic, unidirectional distri- 
bution of streaming cos mic rays, the rates for the reso- 
nant growth of Alfvenic jPohl fc Schlickeiserl [2000D and 
electrostatic (jPohl et al.ll2002f ) turbulence have been de- 
rived using quasilinear theory. Based on an analytical 
treatment, IReville et al.l (|2006) found that also in this 
case nonresonant, purely growing modes may be ex- 
pected to be significantly faster, although the growth 
rate falls off with the temperature of the background 
medium. Application of this mechanism to the ex- 
ternal GRB shocks was phen omenologically studied by 
iMilosavljevic fc Nakarl (|2006[ L who concluded that CR- 
driven turbulence may account for the levels of amplified 
magnetic fields inferred from these sources. We have per- 
formed a series of two-dimensional simulations for this 
setup to explore the relationship between this instabil- 
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ity and that found for drifting cosmic rays, and to de- 
termine the mutual backreaction between the magnetic 
turbulence and the CR beam. The interaction of a cold 
relativistic ion beam is studied in the limit of a mag- 
netized background plas ma, for which the re sults of the 
analytical calculations of iReville et alJ (j2006) apply. 

The simulation setup is described §2 and the results 
of the linear kinetic analysis of the system are presented 
in §3. In §4 the simulation results are presented. The 
differences in the properties of the system between the 
runs representing the CR beam with a constant external 
current and the fully kinetic simulations are discussed in 
§4.1 and §4.2. The detailed properties of the magnetic 
turbulence and the evolution of particle phase-space dis- 
tributions are then presented in §4.3 based on the results 
for mildly relativistic ion beams. We conclude with a 
summary and discussion in §5. 

2. SIMULATION SETUP 

The code used in this study is a 2.5D (2D3V) version 
of the relativistic electromagnetic part icle code TRIS- 
TAN with MPI-ba sed parallelization (jBunemanl fl993fc 
iNiemiec et al.ll2008l ). In the simulations a cold, relativis- 
tic, and monoenergetic cosmic-ray ion beam with Lorentz 
factor 7c r (velocity vcr) and number density Ncr 
streams along a homogeneous magnetic field Bug relative 
to the ambient electron-ion plasma. The ions of the am- 
bient medium have a thermal distribution with number 
density iVj, in thermal equilibrium with the electrons. 
The electron population with density N e — N{ + Ncr 
contains the excess electrons required to provide charge- 
neutrality and drifts with Vd = vcrNcr/ N e with re- 
spect to the background ions, so it provides a return 
current balancing the current carried by the ion beam. 
We have explored the system in the limit of a magnetiz ed 
background plasma, oj <C (see (|Reville et al.ll2006h V 



Specifically, we assumed r y ma x/^i — 0.2, where 
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is the growth rate of the most unstable nonreso- 
nant mode, fli is the ion gyrofrequency, and va = 
[-B 2 / Ho(N e m e + JVjjBi)] 1 / 2 is the plasma Alfven ve- 
locity The relativistic cosmic-ray populations repre- 
sent very dilute ion beams which we study with den- 
sity ratios Ni/NcR = 50 and 125, and Alfven velocities 
va = c/20 and va — c/50, respectively (for which the 
ratio uj pe /rt e = 4.4 and 11.0, respectively). The simula- 
tions have been performed for the case of an ultrarela- 
tivistic beam with "fcR = 300 and a slower beam with 
7cfl = 20. We have also studied the case in which the 
beam is represented by a constant uniform external cur- 
rent, so that the backreaction of the magnetic turbulence 
on the cosmic-ray beam is suppressed. The parameters 
of all simulation runs described here are summarized in 
Table 1. 

To ensure the numerical accuracy of our simulations, 
we use a total of 16 particles per cell, and we apply the 
splitting method for the beam particles. The density ra- 
tio between simulated cosmic-ray and ambient particles 
on the grid is 1/3, and weights are applied to each beam 
particle to match the desired NcR/Ni. The same weights 
are used for the excess electrons 8N e = N e — Ni = Ncr, 



and thus each ion particle can be initialized at the same 
location as the corresponding electron for the identi- 
cally zero initial charge density. The electron skindepth 
A se = c/ujpe = 4A, where ui pe — (A r e e 2 /m e eo) 1 ^ 2 is the 
electron plasma frequency and A is the grid cell size. 
We further assume a reduced ion-electron mass ratio 
rai/m e = 20. This choice allows us to clearly separate 
the plasma and turbulence scales and yet use a compu- 
tational box that can contain several wavelengths of the 
most unstable mode 



X,, 



(2) 



where X S i is the ion skindepth. As the results of our lin- 
ear analysis of §3 sho w, the one-dimension al dispersion 
relations obtained by IReville et al.l (|2006f ) capture the 
physical properties of the system, provided the Lorentz 
factor of the ion beam is very large or a constant ex- 
ternal current is applied. In this study we perform 
simulations for an ultrarelativistic beam and for the 
case of a constant external cosmic-ray current to cross- 
check those results as well as the validity of our ap- 
proach. For these runs we use smaller computational 
grids: (L x ,L y ) = (7AX max , 5.7X max ) for runs A and 
D, and (L x , L y ) = (A.6X max , 3.3X max ) for run B. Larger 
grids with (L x ,L y ) = (10.2X max ,7AX max ) are used to 
investigate the systems with mildly relativistic cosmic- 
ray beams, which have not been studied with PIC sim- 
ulations before. Cosmic-ray beam ions move in the 
— x-direction, antiparallel to the homogeneous magnetic 
field -B||o- Periodic boundary conditions are assumed 
for all boundaries. In all simulations the time step is 
St = 0.1125/tJp e , and the inverse maximum growth rate 
of the nonresonant modes 7"^ = 39758t for runs A-C 
and 9938<5t for runs D and E. 

Our previous work on similar systems as well as ad- 
ditional test runs performed for the current study en- 
sures that our results are not affected by a particu- 
lar choice of simulation parameters, e.g., electron- ion 
mass ratio, number of particles per cell, or the electron 
skindepth. The validity of two-dimensional simulations 
in capturing the e ssenti al physics has been demonstrated 
in INiemiec et ail (|2008l ). 

3. LINEAR ANALYSIS 

The growth rate and wavelength of the most unstable 
purely-growing nonreson ant mode gi ven b y equations [1] 
and GO were obtained by iBelll (|2004j ) and IReville et ail 
( 2006) using linear kinetic analysis in the limit of a cold 
ambient plasma and only for wavevectors fey parallel to 
Biiq. We have numerically calculated the growth rates 

for arbitrary orientation of the wavevector, k, in the zero- 
te mperature l imit (for details of the calculations see, e. 
g., (Bret 2009)). In Figure 1 we show the growth rates in 
the reduced wavevector space {Z\\, Z±), Zi — kiVcR/u pe , 
that is contained in our simulation box, for beams mov- 
ing along a homogeneous magnetic field of strength given 
by the Alfven velocity of va — c/20. In each case the 
dominant unstable mode is the electrostatic Buneman 
mode between background ions and drifting electrons 
(Zii > 2). The growth rate of this mode is about 10 2 
times larger, and its wavelength about 1.25 x 10 3 shorter 
than that of the nonresonant mode (see Eqs. [T]and[2]). 
Our simulations do not fully resolve this mode, because 
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the wavelength of maximum growth corresponds to half 
a cell (Z|| ~ 55) on our computational grid, and we are 
able to see longer wavelength modes only with somewhat 
smaller growth rates. However, the Buneman instability 
is very sensitive to thermal effects and should saturate if 
the thermal velocity of ambient particles becomes compa- 
rable to their relative drift velocity. The initial electron 
thermal velocity in the simulations is thus set to values 
v e .th ^ Vd to ensure the quick saturation and dissipation 
of this unstable mode. Note that such plasma parame- 
ters well reproduce the real conditions in astrophysical 
objects, and the Buneman mode will be relevant only if 
the beam density is high, because then Vd will be high as 
well. 

The nonresonant mode, which we are chiefly interested 
in, is visible at Z±_ < 0.1 and shows a broad peak centered 
at Zu w 0.05, corresponding to the estimate given by 
Eq.[2] However, if the CR beams are treated fully kinet- 
ically (jcr = 300 and 20 - Fig. [lb and Fig. [lb, respec- 
tively), the nonresonant mode is not dominant even in 
the limited wavevector space covered in our simulations. 
In fact, the strongest growth occurs for 0.3 < Z± < 1, 
almost independent of Zu . The very peaked growth at 
Zu w 1 pertains to the Buneman instability between rel- 
ativistic ion beam and ambient electrons. The growth at 
smaller Z\\ represents the filamentation of the ambient 
plasma and the ion beam. The appearance of these fast- 
growing modes modifies the system, and one should ex- 
pect that the properties of the nonresonant mode emerg- 
ing in the nonlinear stage diff er from those predic ted in 
the analytical calculations bv lReville et a l. (2006). This 
is in fact what is observed in our simulations. It should 
be noted, though, that for warm CR beams filamentation 
might be suppressed, and therefore may not play a role 
in the precursor region close to the shock. 

Note that the growth rates as shown in Figure 1 depend 
on the parameters of the system under study. In partic- 
ular, for ultrarelativistic beams (Fig. [Tb>) the growth 
rate for the filamentation modes is much smaller than 
for 7c/; = 20. Hence, for ^cr. ^ 100 we primarily have 
a competition between the nonresonant mode and the 
Buneman modes. If we replace the ion beam by an con- 
stant external current (which somewhat corresponds to 
7 3> 1, Fig. QJ,), then the ambient electrons do not inter- 
act with the ion beam, and a Buneman instability is not 
excited. The evolution of the system is the n artificially 
domi nated by the nonresonant instability l |Ohira et alj 
12001 . 

4. SIMULATION RESULTS 

4.1. Simulations with Constant Cosmic-Ray Current 

The temporal evolution of the energy density in the 
transverse magnetic-field component is shown in Fig- 
ure [2k.. If the backreaction on the cosmic rays is sup- 
pressed, i.e., a constant uniform external current is 
applied, then a purely-growing parallel mode of mag- 
netic turbulence appears in the plasma. Its growth 
rate (7 « 0-Sjmax for the two cases with va = c/20 
and c/50, run A and D, respectively) and wavelength 
(dashed line in Fig. [5b) agree well with those predicted 
by quasi-linear analytical calculations. The mode repre- 
sents a purely magnetic, circularly polarized, and aperi- 
odic transverse wave. The interactions of the magnetic 



turbulence with the plasma are predominantly related to 
the return current carried by the ambient electrons, jret- 
The j re t x 5B± force induces motions and turbulence in 
the background plasma, which in the later stages cause 
the turbulence to turn nearly isotropic and highly nonlin- 
ear. A s in the case of drifting cos mic ravs (INiemiec et akl 
I2008D and nonrelativistic beams (jWinske fc Lerovil 19841) . 
the saturation of the magnetic-field growth proceeds via 
bulk acceleration and occurs when the bulk velocity of 
the background plasma approaches the cosmic-ray ion 
beam speed (see §4.3.2.). Note that current and charge 
balance is still observed if the cosmic-ray current is cho- 
sen constant, and the background plasma is charged by 
adding extra electrons to compensate for the charge of 
the cosmic rays. The plasma thus "knows" the cosmic- 
ray drift speed as that at which the plasma electrons 
no longer stream relative to the plasma ions to carry the 
return current. The nonlinear amplitude of the field per- 
turbations is sli ghtly larger for smaller A lf yen v elocity, in 
agreement with iRiquelme fc Spitkovskvl (|2009f ). How- 
ever, the magnetic-field amplitudes become comparable 
at the end of both runs, and reach SB±/Bu ~ 25, which 
is close to the maximum obtained with M HD simulations 
(|Bellll2004 [20051 iZirakashvili et al.ll2008h and other PIC 
simulations ( Ohira et al.ll2009l ). in which the cosmic rays 
were also represented by a constant current. 

We will now describe the behavior of the system includ- 
ing the response of the relativistic cosmic-ray ion beam. 

4.2. Fully Kinetic Simulations 

If the cosmic rays are treated fully kinctically, the dy- 
namics of the system changes. The interaction of the ion 
beam with the plasma quickly leads to plasma and beam 
filamentation which is modified by a Buneman instability 
between the ion beam and plasma electrons. The Bune- 
man beam-electron interactions produce mainly electro- 
static, slightly oblique turbulence whose wavelength par- 
allel to the direction of the beam is in very good agree- 
ment with the predictions of our linear analysis, which 
gives A = 2tt{vcr/ c)\ se ~ 25A (§3). The mode grows 
very fast, causing density fluctuations in the beam and 
electron plasma. However, in the simulations its am- 
plitude quickly saturates and is subsequently kept at 
a moderate level. These features are in agreement 
with the known properties o f the Buneman modes (see, 
e.g.. lDieckmann e~ al. (2007) for a detailed discussion of 
the nonlinear evolution and saturation mechanism of the 
Buneman instability). The Buneman mode dissipates 
only after filamentation and nonresonant modes have 
strongly backreacted on the ion beam in the nonlinear 
stage. 

4.2.1. Ultrarelativistic Beams 

The properties of the magnetic turbulence depend on 
the Lorentz factor of the beam. For an ultrarelativistic 
beam with jcr — 300 (run B; dotted lines in Fig. 2), the 
filamentation is weak and the parallel nonresonant mode 
appears with the theoretically predicted wavelength. Its 
growth rate is initially 7 ss 0.947 maa; and decreases dur- 
ing the nonlinear evolution. As one can see in Figure 2a, 
the peak amplitude of the magnetic-field perturbations, 
5B±/B\\ a ~ 9, is close to that obtained with constant ex- 
ternal current (run A; solid line) at the onset of the satu- 
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ration of the turbulence growth (t ~ 157"^). It appears 
that in this phase the high beam Lorentz factor provides 
sufficient stiffness to the ion beam that its backreaction 
is suppressed, rendering the system response similar to 
that for a constant external current. The similarity ends 
when the saturation kicks in, though. The subsequent 
dissipation of the turbulence in the run with ^cr = 300 
is much stronger than in the case of a constant external 
current, which places in doubt the accuracy of simula- 
tions that use a constant external current to describe the 
highly nonlinear phases in the evolution of the system. 

4.2.2. Mildly Relativistic Beams 

Results for a system with a mildly relativistic beam 
with jcr = 20, and for va = c/20 (run C) and c/50 
(run E), are presented in Figure 2a with dash-dotted 
and long-dashed lines, respectively. As our linear analy- 
sis of §3 shows, the filamentation modes at perpendic- 
ular wavevectors k± m 1/A se are strong in this case. 
They cause filamentation in the ambient plasma and the 
ion beam, before the nonresonant parallel modes have 
emerged. As one can see in Figure 2, these modes do not 
lead to magnetic-field perturbations of significant ampli- 
tude. Nevertheless, their action on the ambient plasma 
changes its properties, which considerably influences the 
characteristics of the purely-growing parallel modes. The 
nonresonant modes appear in a broad range of wave- 
lengths around X rna x (Fig. 2b), and the growth rate of 
the magnetic-field perturbations is only ~ 0.4r/ max . The 
backreaction of the turbulence on the system further en- 
hances the filamentation in the beam and the plasma, 
and leads to the saturation and dissipation of the mag- 
netic turbulence at a level a few times the homogeneous 
magnetic field strength. The peak amplitudes for the 
two cases with va = c/20 and c/50 are 8Bj_/Bn ~ 4.7 
and 5B±/B\\ ~ 7.5, respectively, showing that instabil- 
ities operating in a less-magnetized medium provide a 
stronger field amplification. It is unclear whether the 
modification of the parallel mode arises specifically from 
filamentation or from any type of perpendicular, small- 
scale density fluctuations, including preexisting turbu- 
lence. We can therefore not reliably predict the behavior 
of a system containing a warm cosmic-ray beam, for a 
example the denser parts of a cosmic-ray precursor to an 
astrophysical shock. 

4.3. Aperiodic Magnetic Turbulence Produced by Mildly 
Relativistic Beams 

4.3.1. Spectral Properties of the Turbulence 

The characteristic features of magnetic turbulence in a 
system containing a mildly relativistic cosmic-ray beam 
are detailed in Figures El [U and [5] for the run with 
-fCR = 20 and va = c/20 (run C). The temporal evo- 
lution of the magnetic and electric field average energy 
densities is shown in Figure [3J Figure [H presents Fourier 
power spectra of the perpendicular magnetic-field com- 
ponent B z for t'jmax — 2, 5, and 8 in two-dimensional 
reduced wavevector space (Zji, Z±). Figure [5] shows snap- 
shots of the time evolution of the electron and cosmic-ray 
ion density, and the structure in the B z magnetic-field 
component. 

The initial filamentation in the ambient plasma grows 
quickly in spatial scale by merging of adjacent filaments, 



which can be clearly seen in E x and E y , and also in 
the B z field components shown in Figure |U Because 
the Buneman instability between the ion beam and the 
plasma electrons is slightly oblique (see Fig. [Tfc), hence 
not purely electrostatic, it is visible in magnetic-field 
Fourier spectra as a feature at Z\\ » 1 (Fig. [4j. The 
corresponding strong short-scale modulations in the den- 
sities of ambient electrons and the ion beam can be seen 
in Figure [5^,-b and Figure [5ji-e. The nonresonant par- 
allel modes of magnetic turbulence emerge in a medium 
already strongly modified by filamentation. They ap- 
pear in a range 0.02 < Z» < 0.1 around the theoretically 
predicted Zn{X max ) ~ 0.08 and quickly grow in wave- 
length (see Figs. |3Jd-c,[5J:, and[Sf). The influence of the 
nonresonant modes is stronger on the filamentation in 
the slowly drifting ambient plasma than that in the rel- 
ativistic ion beam. In essence, ambient plasma filaments 
become vertically tilted (Fig. [5U), which leads to even 
stronger plasma filamentation. The lack of spatial corre- 
lation between filaments in the ambient plasma and the 
beam results in a local charge imbalance and the build- 
up of charge-separation electric fields, which, together 
with electric fields induced by the Buneman instability, 
dominate the turbulent electromagnetic energy content 
of the system in the initial stage (Fig. [3]). During the 
nonlinear stage (t > 87"^) the enhanced filamentation 
leads to the generation of stronger turbulence in the 
B z component of the magnetic field with perpendicu- 
lar wavevectors, k y , which disrupts the structure of the 
parallel magnetic modes. This interaction between fila- 
mentary and nonresonant modes is visible in Figure [Sf 
and in the Fourier spectrum in Figure [4);. 

As one can see in Figures [5^-i, the strongly amplified 
magnetic field starts to backreact on the cosmic-ray beam 
in the later stage of the system evolution. Cosmic-ray 
filaments become tilted and eventually disrupted. This is 
accompanied by turbulent ambient plasma motions and 
results in highly nonlinear and nearly isotropic magnetic 
turbulence. The characteristics of the turbulence in its 
post-saturated state are thus similar to those observed in 
simulations of nonrelativistically drifting hot cosmic-rays 
in the precursor to SNR shocks. 

4.3.2. Particle Phase-Space Distributions 

The effects of the backreaction of the magnetic turbu- 
lence on the particles are presented in Figures [7J and 
[SI The average (bulk) velocities of all particle species 
converge in the nonlinear stage, when the magnetic-field 
growth saturates (at t ~ 157"*,,. for runs A and C in 
Fig. [6]). While the relative drift between the plasma 
and the cosmic-ray beam disappears in all our simula- 
tions, the mechanism by which that is achieved differs 
between runs which suppress the cosmic-ray backreac- 
tion and fully kinetic runs. 

In the fully kinetic simulations (solid line in Fig. 0, 
the cosmic-ray beam slows down considerably, and the 
ambient plasma accelerates up to ~ —0.2c. This beh av- 
ior is similar to the results of Winske fc Lerovl (|1984D for 
nonrelativistic dilute ion beams interacting with ambient 
plasma via nonresonant modes, which showed that the 
energy of the decelerating beam is transferred in approx- 
imately equal parts to the ambient ions and the magnetic 
field, in accordance to the predictions of quasi-linear the- 
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ory. The simulation results of iWinske fc Lerovl (|1984f ) 
were obtained with a one-dimensional hybrid model that 
treats electrons as a massless fluid, with which one can- 
not observe filamentation modes, and the associated elec- 
tron heating is artificially suppressed. As shown in Fig- 
ure which presents the temporal evolution of energy 
densities in particles and fields for run C, in our simula- 
tions the initial plasma filamentation is accompanied by 
electron heating at the expense of the beam. However, 
the heating of the electrons saturates at t w fjmax' when 
the nonrcsonant modes start to emerge. During the sub- 
sequent evolution, beam energy is transferred with ap- 
proximately the same rate into the magnetic field and 
ambient ions, while the electrons experience only moder- 
ate further heating. This process of the energy transfer 
saturates when the turbulent magnetic field reaches its 
maximum energy density and starts to dissipate. The 
nonlinear evolution of nonresonant modes in a system 
containing relativistic ion beams thus proceeds in quali- 
tatively the same way as for nonrelativistic beams. 

If the cosmic-ray ion beam is represented by a con- 
stant external current, the energy conservation between 
the beam to the ambient medium is violated. As shown 
in Figure [6] (dashed lines), the saturation of the magnetic 
field growth still comes about by the disappearance of the 
ion beam-ambient plasma relative motion, at which the 
return current is provided without a drift of the plasma 
electrons relative to the plasma ions, but now the ambi- 
ent ions and electrons must assume the constant cosmic- 
ray beam bulk velocity. This implies that energy is con- 
tinuously pumped into the system, and therefore energy 
conservation becomes severely violated in the nonlinear 
stage. Thus the validity of simulations that assume a 
constant cosmic-ray current is limited to the early phases 
in the evolution of the system. 

Figure [8] shows the phase-space distributions of the 
cosmic-ray beam and the ambient ions at t^ max — 0, 7, 
14, and 19 for run C. The early stage of the system evolu- 
tion (t^max ^ 7) is dominated by the Buneman instabil- 
ity modes between the ion beam and ambient electrons. 
The electrostatic fields associated with this mode heat 
up electrons (Fig. [7|) and significantly stretch the beam 
ions distribution along the beam propagation direction. 
At the same time the cosmic-ray beam is heated in the 
transverse direction due to filamentation modes. The 
ambient ions remain unaffected by the Buneman mode 1 
and become only moderately heated in this stage. The 
stretching of the beam ions distribution gradually sat- 
urates at t^ max ~ 7, by which time the nonresonant 
modes have set in and started to strongly backreact on 
the system. 

During the subsequent evolution the beam momentum 
becomes quickly randomized in direction. This random- 
ization is the combined effect of the pinching of ion-beam 
filaments and pitch-angle scattering of the beam parti- 
cles. The latter process becomes more important in the 
highly nonlinear phase (tj max > 10; compare Fig. [SJi), 
during which the filaments start to get disrupted. At 
the same time the ion beam slows down in bulk, and 

1 The phase velocity of the Buneman wave mode between cosmic- 
ray ion beam and electrons is ~ v cr in the ambient ions rest frame. 
Thus the associated electrostatic fields are seen by the ambient ions 
as a high-frequency oscillations. 



by tj max ~ 19 the evolution saturates when the ion 
beam particles have been efficiently pitch-angle scattered 
around a mean (bulk) momentum of ~ —3.8miC. 

The randomization of beam momentum through pitch- 
angle scattering was previously reported for nonrela- 
tivistic beams in conditions allowing for an efficient 
magnetic field amplifica tion through nonresonant modes 
(IWinske fc Lerovl Il98l . Here we have demonstrated 
that these modes can also provide efficient scattering for 
relativistic beams. 

We can estimate the scattering mean free path from 
the time evolution of the phase-space distribution of the 
ion beam, a few snapshots of which are shown in Figure[5J 
Between 1 7 max = 10 and 1 7, nax = 14 the scattering mean 
free path in simulation run C is 

A mfp « 5000A. (3) 

At the same time, the rms amplitude of the turbulent 
magnetic field increases by more than 250%, from about 
B||o to 3.5-B||o- Using the mean of the two numbers, we 
obtain for the Bohm mean free path 

A B ohm« 3000A (4) 

Given the uncertainty in the estimate arising from the 
substantial variation in the magnetic-field amplitude, 
about a factor 2, we conclude that the observed scat- 
tering mean free path, and therefore the spatial diffusion 
coefficient, for mildly relativistic beams are entirely com- 
patible with Bohm diffusion. 

The estimate of the scattering mean free path can also 
be made for ultrarelativistic beams based on run B. How- 
ever, by the time our simulation ends the CR beam is 
only partially pitch-angle scattered up to an angle ~ 7r/6. 
Nevertheless, a rough estimate shows that A m f p is again 
within the factor of a few comparable to Aeohm- 

5. DISCUSSION AND CONCLUSIONS 

We have studied the interaction of a cold, relativistic 
ion beam penetrating a cold plasma composed of elec- 
trons and ions. We have presented 2.5D PIC simula- 
tions, complemented with a linear analysis of the dis- 
persion relation for linear waves with arbitrary orienta- 
tion of k, for parameters that permit the g rowth of non- 
reson ant, purely-magnetic parallel modes (jReville et al.1 
2006). Our research is relevant for the understanding of 
the structure of, and particle acceleration at, shocks in 
SNR, GRB, and AGN, for which radiation modeling sug- 
gests that the magnetic field near the shock is strongly 
amplified. 

We observe a close competition of the nonresonant 
mode with the filamentation instability and Buneman 
modes, which is also evident in the linear dispersion rela- 
tion. The specific choice of parameters determines which 
of the three modes of instability dominates. In some 
cases filamentation is initially important and modifies 
the later evolution of the parallel nonresonant mode. In 
all cases we find that a representation of the ion beam by 
a constant current, as is routinely done in MHD studies, 
is suboptimal, because it suppresses part of the nonlin- 
ear response of the system, delays the saturation pro- 
cesses, and leads to a significant overestimate of the final 
magnetic-field amplitude. 

As i n the case of drifting cosmic rays (Niemiec et al.l 
120081: IStroman et a l. 2009) and nonrelativistic beams 
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(jWinske fc Lerovlll984f) . the saturation of the magnetic- 
field growth proceeds via bulk acceleration. For mildly- 
and ultra-relativistic beams, the instability saturates 
at field amplitudes a few times larger than the ho- 
mogeneous magnetic field. These results match our re- 
cent studies of nonrelativistically drifting cosmic-rays up- 
stream of SNR shocks which also indicated only a mod- 
erate magnetic-field amplification by nonresonant insta- 
bilities. 

We have demonstrated that the magnetic field ampli- 
fied via nonresonant interactions between the CR beam 
and the plasma can efficiently scatter cosmic rays even for 
moderate field amplification levels. The scattering mean 
free path is compatible with Bohm diffusion. Sub-Bohm 
diffusion was observed in Monte-Carlo simulations of par- 
ticle transport in the nonlinear turbulent magnetic field 
genera ted in the nonresonant instability bv lReville et al.l 
(2008). In that work, parallel and perpendicular diffu- 
sion coefficients were calculated by probing the spatial 
displacement of test particles in a static snapshot of the 
amplified magnetic field that resulted from MHD simu- 
lations of the instability. Here we estimate the isotropic 
spatial diffusion coefficient by probing the evolution of 
the angular distribution of particles in the self-excited, 
non-stationary (growing) turbulence whose typical wave- 
length is at least a factor of a few smaller than the gyro- 
radii of cosmic rays. 

In the application to nonrelativistic shocks in SNRs, 
strong (SB /Bo 3> I) quasi- isotropic magnetic turbulence 
would be compressed by the shock, thus turning into 
quasi-two-dimensional turbulence in the downstream re- 
gion. Radio polarimetry suggests that the magnetic field 
immediately behind the shock is preferentially oriented 
along the shock normal (jStroman fc Pohill2009f ). which is 
at odds with the above expectation, if the tur bulent field 
is not quickly damped to an amplitude < B» (jPohl et al.l 
l2005h . 
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ing sufficient to randomize CR momenta. However, it 
is not clear how the strong quasi-isotropic magnetic tur- 
bulence in the upstream region influences the parti cle 
acceleration at the shock (but see lCouch et al.l (|2008l) ). 

Finally, our simulations show that the saturation of 
instabilities operating upstream may limit the magnetic 
amplitude to moderate levels. If very strong magnetic 
field is required by radiation modeling, it may there- 
fore be generated at the shock itself or in the immediate 
downstream region. 
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Fig. 1. — Linear growth rate 7 in units of the electron plasma frequency, LOp e , as function of the flow-aligned and perpendicular (fey and 

fcj_) wavevectors for the parameters used in simulations A, B, and C (see Table 1). For a given fc, the growth rate of the most unstable 
mode is plotted. The figures show the modes whose wavelengths are well contained in the simulation box. 
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Fig. 2. — Temporal evolution of the energy density in the transverse magnetic field component, normalized to the homogeneous field 
strength (a). Time is dimcnsionless in units of the theoretically predicted inverse growth rate of the nonresonant mode, Jmax- Simulations 
with constant cosmic-ray current are presented with the solid line for a plasma Alfven velocity va = c/20 (run A), and the dashed line for 
va = c/50 (run D). The effects of cosmic-ray backreaction are shown for the case of va = c/20 and the beam Lorentz factor -ycR = 300 (run 
B; dotted) and ~/cr = 20 (run C; dashdottcd), and for va = c/50 and "fCR = 20 (run E; long dashed). Fourier spectra of the perpendicular 
magnetic-field component, B z , in wavelengths along the beam direction for the case with va = c/20 (runs A-C) (b). The parallel mode 
is expected at the wavelength X m ax- The dashed line shows results for a constant cosmic-ray current (run A), and the dotted line for 
a faster beam with "fCR = 300 (run B; spectrum shifted for presentation clarity) at the initial stage of the turbulence growth, ty max = 
7 and 9, respectively. Note the single peak structure centered on \ m ax- The solid lines show the temporal evolution of the spectra for 
■fCR = 20 (run C) during the linear stage (t-y ma x = 6, 7, 9, 11, from bottom to top). The magnetic turbulence appears in a broad range 
of wavelengths. 



TABLE 1 

Simulation Parameters and Results 



Run 


Grid 


^rnax 


Ni/NcR 


ICR 


VA 


'j/ymax 


6B™ a */B ll0 




i^max ) 


{.Imax) 






(c) 






A 


7.4x5.7 


32.7 


50 


oo 


1/20 


0.83 


24.0 


B 


4.6x3.3 


21.0 


50 


300 


1/20 


0.94 


8.9 


C 


10.2x7.4 


31.2 


50 


20 


1/20 


0.43 


4.7 


D 


7.4x5.7 


30.2 


125 


oo 


1/50 


0.79 


24.0 


E 


10.2x7.4 


26.7 


125 


20 


1/50 


0.4 


7.5 



Note. — Parameters and selected results of the simulation runs described 
in this paper. Listed are: the system size in units of A^^,, the run duration 
t max in units of 7^„, the density ratio of ambient and beam ions, beam 
Lorentz factor, the Alfven velocity in units of the speed of light c, the measured 
growth rate 7 of the nonresonant modes in units of 7 ma;E , and the maximum 
amplitude of the perpendicular magnetic field perturbations 8B™ ax relative 
to the homogeneous magnetic field. The other parameters assumed in the 
simulations are: the electron skindepth A se = 4A, the ion skindepth A s ^ ~ 
18 A, \max = 562 A, ■y m ax/^li = 0.2, and the ion-electron mass ratio m,i/m e = 
20. The runs assuming a constant cosmic-ray current are marked with -fCR = 
oo. The growth rate 7 is provided for time periods in which nonresonant 
parallel modes show up in the plasma. The amplitudes of the field turbulence 
for runs with constant cosmic-ray current (A and D) represent the unsaturated 
values at the end of these simulations. 
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Fig. 3. — Temporal evolution of the average energy density in electromagnetic fields, normalized to the initial total kinetic energy in the 
system, Uko, f° r the run with 7ci? = 20 and Va = c/20 (run C; compare Fig. [2p). Here B^ and indicate the magnetic and electric 
energy densities in the transverse components, i.e., {By + B^) / (2/j.o) and correspondingly for the electric field. The quantities Bp <5-Bjj , 
and E? are defined analogously, but SB?, refers to the turbulent component only. 




Fig. 4. — Fourier power spectra [log 10 (dP / d\)] of the perpendicular magnetic-field component B z in the run with ■yen = 20 and Vj\ = c/20 
(run C) for t = 2,5, and 87^0^ in two-dimensional reduced wavevector space (Z^,Z±). The wavevector range is chosen to facilitate a 
direct comparison with analytically predicted dispersion relation shown in Fig. [l] 
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Fig. 5. — Snapshots of the time evolution of the ambient electron density (a, d, g), cosmic-ray ion density (b, e, h), and the B z field 
amplitude (c, f, i) for run C at t"f max = 6 (top), 9 (middle), and 14 (bottom). Electron and ion beam densities are normalized to their 
respective initial values. The magnetic field B z is normalized to the amplitude of the homogeneous field, Bu . The spatial scales are 
provided in units of X m ax and only a small portion of the simulation box ((Z/ x , L y ) = (10.2A maI , 7 AX ma x)) is shown. 
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Fig. 6. — Temporal evolution of the bulk (average) velocity of all particle species for run A (dashed lines) and run C (solid lines). 
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Fig. 7. — Temporal evolution of the average energy density in cosmic-ray beam ions, Ucr, ambient electrons, U e , and ions, Ui, and 
the volume-averaged electromagnetic energy density, Uem = ((-^x + By + B^)/2/j,o + (E% + Ey + E^)/2), for run C. All quantities are 
normalized to the initial total kinetic energy in the system, Uko- 
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Fig. 8. — Phase-space distributions of the ions in the cosmic-ray beam (gray) and the ambient plasma (black) at tjmax = 0, 7, 14, and 
19 for the run with -fen = 20 and va = c/20 (run C). 



